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We investigate dynamics of activated escape in periodically modulated systems. The trajectories 
followed in escape form diffusion broadened tubes, which are periodically repeated in time. We 
show that these tubes can be directly observed and find their shape. Quantitatively, the tubes are 
characterized by the distribution of trajectories that, after escape, pass through a given point in 
phase space for a given modulation phase. This distribution may display several peaks separated by 
the modulation period. Analytical results agree with the results of simulations of a model Brownian 
particle in a modulated potential. 
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I. INTRODUCTION 

The problem of noise-activated escape in periodically 
modulated systems is of interest to many areas of physics 
and applications, from Josephson junctions 0,0,13 and 
nano- and micromechanical systems 0, to epidemics 
0. The theory of activated escape should answer two 
closely related questions: what is the escape rate and how 
does the system move during escape? Fluctuational tra- 
jectories leading to escape or interstate switching are of 
significant interest and have been extensively studied in 
recent years using various numerical techniques, see Refs. 
IE M and papers cited therein. Even though escape 
is a random event, the probabilities of following differ- 
ent paths are strongly different. Therefore most likely 
the system follows a certain pathway, i.e., its trajectory 
is close to the most probable escape path (MPEP), see 
Ref. E and references therein. 

The distribution of escape trajectories can be conve- 
niently characterized by the prehistory probability distri- 
bution (PPD), a quantity accessible to direct experimen- 
tal measurements. It is obtained by recording trajecto- 
ries of the system that lead to escape and superimposing 
the trajectories that, after the system has escaped, pass 
through a small vicinity of a point qj (for a certain mod- 
ulation phase, in the case of a periodically modulated 
system). The point g/ is chosen in the phase space suffi- 
ciently far behind the boundary of the basin of attraction 
to the initially occupied metastable state. Formally, the 
PPD ph(q,t\qf,tf) is the probability density for the sys- 
tem to have passed through a point q at an instant t 
provided it had been fluctuating about the metastable 
state for a long time and passed qf at a later time tf, 

t f >t\u&. 

The PPD ph(q,t\qf,tf) should peak at q lying on the 
optimal fluctuational path that leads to qf. Therefore it 
"maps out" optimal paths. For stationary systems this 
has been directly confirmed by extensive simulations fl^ . 
E,E] and also in laser experiments In such systems 
escape can occur at any time, with equal probability, 
therefore the tube of paths around the optimal escape 
path is broad. 

In the present paper we study escape pathways in pe- 



riodically modulated systems. Modulation synchronizes 
escape events. This can be easily understood for escape 
from a slowly modulated potential well. Here, escape 
is most likely to occur when the instantaneous poten- 
tial barrier AU(t) is at its lowest, once per period, cf. 
Fig-da). If the modulation amplitude of AU(t) signifi- 
cantly exceeds the noise intensity D (D — ksT for ther- 
mal noise) , escape events are synchronized very strongly. 
Strong synchronization persists even where the modula- 
tion is not slow and the adiabatic picture in which the 
escape rate is determined by the instantaneous barrier 
height does not apply E H3, E E E] As a result of 
escape synchronization, there is one MPEP per period. 

In turn, as we show, the prehistory probability distri- 
bution displays a sharp narrow peak as a function of q 
for a given t. This peak lies on the MPEP. Respectively, 
in (q, t) space the PPD displays a narrow ridge centered 
at the MPEP, cf. Fig. [3 Moreover, the PPD may have 
several narrow peaks [ridges in (q, t) space] . Their width 
gives the width of the distribution of escape trajecto- 
ries. It is determined by the typical diffusion length 
Id = (2-Dir) 1 / 2 , where t x is the relaxation time of the 
system. This is qualitatively different from the shape of 
the PPD in stationary systems El- Of course, in modu- 
lated systems along with the final point qf through which 
the system passes one should fix the modulation phase 
when the passage happens. It is given by the passage 
time i/(mod Tp ), where Tp is the modulation period. 

The occurrence of a narrow peak of the escape trajecto- 
ries distribution can be understood from the qualitative 
picture of motion in escape. This picture is sketched in 
Fig.^for a system with one dynamical variable q. Escape 
from a static potential well corresponds to going over the 
barrier top qb from the vicinity of the potential minimum 
q a , see Fig.QJa). Similarly, a modulated system escapes 
when it goes over the periodic basin boundary qb{t) from 
a periodic metastable state q a (t), see Fig.2Jb) (here and 
below we assume that the noise correlation time is small; 
see Ref. |^| for a more general case) . Let us suppose that 
the escaped particle is found at time tf at a point qf 
sufficiently far behind qb(t). A typical trajectory to this 
point displays four distinct sections with different types 
of motion shown schematically by letters A through D. 
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FIG. 1: (Color online) Activated motion leading to detection 
of an escaped particle at point qf at time tf (schematically). 
Panels (a) and (b) illustrate escape from a static and a peri- 
odically modulated potential well, respectively. In the latter 
case q a {t) and qb(t) are the periodic stable state and the basin 
boundary. The trajectories ?opt(t) ( n — 1, . . . , — 2) are the 
periodically repeated most probable escape paths. The four 
major stages of motion A, B, C, and D on the way to q; are 
discussed in the text. 



We discuss the motion backward in time from t / . Im- 
mediately adjacent to qf is section A of the trajectory in 
Fig. ^ Here, for small noise intensity the system moves 
close to the noise-free trajectory from the vicinity of tft(i) 
to qf. In the case of a static potential this corresponds to 
sliding down the potential slope from qb to qf. Section B 
corresponds to diffusion in a region of width ~ I d around 
qb (t). In this region noise- free motion with respect to 
qb(t) is slow (in the case of a static potential the poten- 
tial is locally flat at qb). Section C corresponds to motion 
from the attractor to the basin boundary. The respective 
trajectory is close to the most probable escape path. This 
motion is a result of the large fluctuation that has led to 
escape detected at qf. For strongly synchronized escape, 
there is one MPEP per modulation period, as shown in 



Fig-Gib). The MPEP approaches qb(t) asymptotically as 
t — > oo. If the system was observed behind qb{tf) at time 
tf, it has most likely arrived to the vicinity of qb (t) along 
the MPEP that approached tfa(i) before tf, but not too 
much in advance, as shown in Fig. ^b). Finally in the 
region D well before t f the system was fluctuating about 
the attractor. 

The above picture suggests that the PPD will peak on 
the trajectory singled out in Fig. ^ It also explains why 
the PPD peak along q axis should be narrow: in contrast 
to stationary systems, which can arrive to the vicinity of 
qb at any time, periodically modulated systems approach 
qb{t) only once per period, and the tube of escape trajec- 
tories is narrow. We note that, by changing qf,tf, we can 
switch between neighboring MPEP's, and therefore there 
is a possibility for the PPD to have two and potentially 
even more peaks inside the attraction basin. 

As mentioned above, we are interested in the regime 
where escape events are strongly synchronized by the 
modulation. It means that the probability density to 
find the escaped system at a point qf behind the basin 
boundary, cf. Fig. ^ displays sharp peaks as a func- 
tion of the observation time tf. These peaks are peri- 
odically repeated in time. They were studied in Ref. I2JJ. 
We will analyze the PPD for the time tf close to the 
maximum of the escape probability peak. This choice 
is justified, because the corresponding PPD character- 
izes the most probable escape trajectories. Otherwise 
the PPD would be formed by trajectories with exponen- 
tially smaller probabilities that are very rarely followed 
in escape. 

Interestingly, the condition that the tubes of escape 
trajectories be narrow in time and space requires that 
the modulation frequency cop — 2tt/tf lie within a range 
limited both from below and from above. To understand 
the lower limit we note that, as mentioned above, for slow 
modulation escape occurs every period around the time 
t m + krp where the height of the instantaneous potential 
barrier AU(t) is at its minimum (k = 0, ±1, ...). The 
typical width of the time window for escape At is given 
by the condition |AZ7(f ro ± At) - AU(t m )\ < D, which 
leads to At = [D/AU(t m )}^ 2 ~ ^[D/AU^)] 1 / 2 . 
If At exceeds the typical duration of motion in escape 
t r \n[AU(t m )/ D], the PPD has the same shape as if the 
system were escaping out of a stationary potential well of 
height AU(t m ). In this case inside the attraction basin 
the ridge of the PPD in (q, t) space is broad and asym- 
metric, its width along q axis is independent of the noise 
intensity |l5|. 

In the opposite limit of high-frequency modulation, 
LUpt T 3> 1, escape of an overdamped system is not syn- 
chronized. The dynamics is characterized by the coordi- 
nates averaged over modulation period. The PPD as a 
function of such coordinates is described by the theory 
for stationary systems [l^ , and the PPD peak inside the 
attraction basin is broad. 

In the intermediate range of frequencies not only are 
the PPD peaks narrow but, as mentioned above, the PPD 
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FIG. 2: (Color online) The prehistory probability density 
(PPD) Ph(q, t\q/,tf) and its contour plot for a noise-driven 
overdamped system with equation of motion q = q 2 — 0.25 + 
AcosuJFt + f(t), where /(t) is white noise of intensity D. 
The parameters are A = 0.7, ujf = 2, D = 0.01, g/ = 0.8, 
tf — [tf /2) (mod Tf), where tf = 2tt/ujf is the modulation 
period. The shadowing (color code on line) corresponds to 
the 4 regions of the height of the distribution separated by 
the values ph = 0.5, 2, 7. 

as a function of q for given t (and also of t for given 
q) may display several peaks inside the attraction basin. 
This happens because the motion of the system near the 
basin boundary is slow. Therefore if the system is ob- 
served behind the boundary qb(t) at a given time tf, it 
could have arrived to the boundary along one of a few 
periodically repeated optimal escape paths, fluctuated 
about qb(t) for some time, and then made a transition 
to {qj,tf) over time ~ t T . 

The multiple-peak structure of the PPD is a specific 
feature of periodically modulated systems far from equi- 
librium, i.e. away from both the adiabatic limit of slow 
modulation and the limit of fast modulation. It is illus- 
trated in Fig. |21 for a model system. Two ridges of the 
PPD inside the attraction basin are clearly resolved in 
this figure. Their shape is well described by the asymp- 
totic theory developed in this paper. 

In Sec. II below we discuss the dynamics of a period- 
ically modulated system and the equations for the most 
probable escape path. In Sec. Ill we obtain a general 
expression for the PPD near the basin boundary and re- 
late it to the distribution that describes the quasiperiodic 
current away from the attraction basin, which gives the 
escape rate. In Sec. IV it is shown that the PPD is sim- 
plified in the adiabatic limit where the modulation period 
is small compared to the relaxation time. Two types of 
behavior may be displayed near qb(t) in this case, de- 



pending on the ratio of the two small parameters that 
characterize the dynamics. The analysis is extended to 
the case of nonadiabatic driving in Sec. V. In Sec. VI 
we obtain the central result of the paper, the PPD in- 
side the attraction basin. We show that it has the shape 
of a sum of diffusion-broadened Gaussian peaks centered 
at the periodically repeated most probable escape paths. 
In Sec. VII we describe the results of simulations of a 
model system and compare them to the analytical re- 
sults. Sec. VIII provides a brief summary of the results. 

II. ESCAPE OF A PERIODICALLY 
MODULATED SYSTEM 

We investigate the prehistory probability distribution 
for an overdamped system characterized by one dynami- 
cal variable q. The system is in a periodically modulated 
potential U(q, t) — U (q, t-\-Tp). Fluctuations are induced 
by an external noise f(t). The motion is described by the 
Langevin equation 

q = K(q,t) + f(t). (1) 

Here, K(q, t) = d q U (q, t) is the regular periodic force. We 
consider the simplest case where f(t) is zero- mean white 
Gaussian noise with correlator (f(t)f(t')) = 2DS(t — t 1 ). 

In the absence of noise the system ^ has a periodic 
metastable state (attractor) q a {t) = q a (t+Tp) and a peri- 
odic boundary qb(t) = qb{t+Tp) of the basin of attraction 
to q a {t), see Fig.^b). The states q a {t) and are, re- 
spectively, the stable and unstable periodic solutions of 
the equation q = K(q,t). For concreteness we assume 
that q a (t) < q b (t). 

The PPD Ph(q,t\q/,tf) as defined in the Introduction 
is the conditional probability density (with respect to q) 
of passing through a point (q, t) in the coordinate-time 
space on the way from the attractor to a point (qf,tf). 
It has the form 11] 

„ r„ ,|„ t \ p(qf,tf\q,t)p(q,t\q in ,t in ) 

Ph{q,t\qf,tf) = / — —, —: , (2) 

PW,tf\Qixi,tin) 

where p(qi,t\\q2,t2) is the probability density (with re- 
spect to qi) of a transition from (92,^2) to (<fi,ti), with 
ti > t%. 

We assume that the noise intensity D is small. Then 
the period-average escape rate is W oc exp(— R/D), 
where R is the activation energy of escape, R ^> D 

[3 IH 13 E The condition on D is that W' 1 

largely exceeds the relaxation time of the system t T and 
the modulation period rp. Eq. gives the distribution 
of escape paths in a broad time range 

W' 1 > tf - t ia >t - t ia > tr, TF) (3) 

where the population of the attraction basin practically 
does not change. For t — t- m t T the initial state q- ln , 
which is close to the attractor (fa (tin), gets forgotten. 
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Then the right-hand side of Eq. J5J) becomes indepen- 
dent of ij n . 

The functions p(q, t\q- ln , t- m ) — p(q,t) and 
p(q f ,t f \q ia ,t in ) = p(qf,tf) give the time-periodic 
probability density to find the system in states q and qf, 
respectively. To study escape pathways we calculate the 
PPD for qf outside the attraction basin, qj > qb(tf). 

The transition probability density p(q,t\q' ,t') is a so- 
lution of the Fokker-Planck equation (FPE) 



d t p = -d q [K(q,t)p]+Dd 2 q p, 



(4) 



with the initial condition p(q,t'\q' ,t') = S(q — q'). Even 
for a ID system, this equation does not have a known 
explicit solution except for the case of K linear in q. To 
analyze the PPD we will have to find approximate solu- 
tions in different regions and match them. 

In the following subsections we discuss the dynamics 
of the system prior and during escape. 



A. Dynamics near the periodic states 

We start with the dynamics close to the periodic states 
qi(t), i = a,b. In the absence of noise it is described by 
the linear equations Sq = pi(t)Sq with Sq — q — <ft(t) and 

=Pt{t + r F ) = [d q K] qi{t) , i = a,b. (5) 
Time evolution of the deviation Sq from qi (t) is given by 
Sq(t) = 5q(t')Ki(t,t'), (6) 



Kj(t, t') = exp 



dTpi(r) 



a, b. 



Because of the periodicity of Pi{t), the Floquet multipli- 



Mi = Ki(t + tjf, t) — exp (imtf) 



a, b, 



(7) 



are independent of time, with M a < 1 and Mj > 1. In 
Eq. pi (i = a,b) is the period-average value of Pi(t), 
with p a < 0, pb > 0. The relaxation time of the system 
can be chosen as t T = p^ 1 ~ | | ~ 1 - 

Weak noise leads to fluctuations about the attractor, 
which have Gaussian distribution near the maximum in 
the time range Q. From Eq. Q this distribution has 
the form 



p a (q,t) =G(q- q a (t);a a (t)) 



(8) 



with 

G(x; a) = (27rL>cr 2 )- 1 / 2 exp(-a; 2 /2Da 2 ). (9) 

The variance cr 2 (t) is periodic in time. It is given by the 
expression 

cr?(t) = 2|Mr 2 -l|" 1 I " dt X Kj 2 {t + tx,t) (10) 
Jo 

with i = a. 



B. Most probable escape paths 

Along with small fluctuations, there also happen occa- 
sional large deviations from q a (t), including escape from 
the attraction basin. In a broad parameter range escape 
events are very strongly synchronized with the modula- 
tion, see Ref. |2jJ and papers cited therein. Because large 
fluctuations have small probabilities, in escape the sys- 
tem closely follows the trajectory that is least improbable 
among all possible escape trajectories. As mentioned in 
the Introduction, it is usually called the most probable 
escape path, q opt (t). Inside the attraction basin, escape 
trajectories lie within periodically repeated narrow tubes. 
The tubes have a width ~ l]j and are centered at the 
MPEP's gg(t) = g opt (f + my), n = 0,±1, . . .. 

The MPEP's provide a solution to the variational prob- 
lem of maximizing the probability of a fluctuation in 
which the system moves from q a (t) to <#,(*). This problem 
can be mapped onto the problem of dynamics of an aux- 
iliary Hamiltonian system. The latter is described by the 
Wentzell-Freidlin Hamiltonian H (p, q; t) = p 2 + pK(q, t) 
|'22j . Its equations of motion have the form 



q = K{q,t) + 2p, p = ~pd q K. 



(11) 



The MPEP's correspond to the minimal action hetero- 
clinic Hamiltonian trajectories (q op t(t),Po P t(t)) [Til l23|. 
They start for t — > — oo at the periodic hyperbolic state 
{la{t),P = 0) of the auxiliary system and for t — > oo ap- 
proach its another periodic hyperbolic state, (qb(t),p — 
0). 

Well outside the diffusion regions around the periodic 
states q a ,b(t) the motion along the MPEP is fast. It 
is seen from Eq. <|11|) that the system moves between 
these regions over time t r . Close to qb(t) the system 
is slowed down, and in the region \q — qb(t)\ < Id the 
motion is dominated by diffusion. The duration of stay- 
ing in the vicinity of qb(t) can be obtained by lineariz- 
ing the equation of motion near qb and is g iven by 
the Suzuki time t$ ~ p^ 1 ln(Ag//u) p^ . |25|| . where 
Aq = min f \qb{t) — q a {t)\ is the typical distance between 
the periodic states. 

Periodically modulated systems are advantageous as 
they allow one to observe, via the prehistory distribution, 
both the fast motion along the MPEP and the slow mo- 
tion near the unstable state. As we show below, the peak 
of the PPD does not display broadening due to diffusion 
near qb, as does the PPD in the absence of modulation 
151. 



III. PREHISTORY PROBABILITY 
DISTRIBUTION NEAR THE BASIN BOUNDARY 

We will calculate the PPD in the regime of strong syn- 
chronization of escape 0, |2(j ■ In this regime the proba- 
bility distribution p(q, t) of finding a particle behind the 
basin boundary has the form of sharp periodic pulses as 
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a function of t. It is most interesting to find the PPD for 
a final point (qf,tf) on the (q, t) plane near the center of 
such a pulse. 

We will assume that the point (q/,tf) is sufficiently far 
from the diffusion-dominated layer (region B in Fig. ^) 
around the basin boundary, so that the distance to the 
boundary is Qf — qf — qb{tf) ^ Id- This condition is 
usually realized in experiments, where the position of a 
particle detector is chosen so as to ensure that the de- 
tected particles have practically no chance to return to 
the attraction basin. At the same time we assume for 
convenience that the final point is still in the harmonic 
region, Qf <C Aq, in which case the motion between qb 
and qf can be described by linearizing equation Q in 
Q = q-qb(t). 

We will start the analysis of the PPD ph(q,t\qf,tf) 
with the case where not only qf but also the point (q, t) 
through which the trajectory of interest has passed is also 
in the harmonic region around the basin boundary qb(t). 



A. Transition probability density 

As seen from Eq. finding the PPD requires cal- 
culating the transition probability density p(qf ,tf\q,t). 
The Fokker-Planck equation for p{qf,tf\q 1 t) can be lin- 
earized near qb(t), 

d tf p = -fi b (tf)d Qf (Q f p)+Dd 2 Qf p, (12) 

where p = p(qf,t f \q,t) and Qf — q f - qb(tf). 

The solution of Eq. 1)12(1 can be sought in the form of 
a Gaussian distribution 



p = G(Qf-Q(tf)-a(tf) 



(13) 



with G(x;a) given by Eq. Q. Eq. (fT2"|l will be satisfied 
provided the hitherto unknown functions a, Q obey the 
equations 

^=2 M6 (t / )a 2 +2, ^- = p b (tf)Q. (14) 

The initial conditions for these equations follow from the 
condition p(qt,t\q,t) = S(qf — q). They have the form 
a 2 (t) = and Q(t) = Q = q — qb(t). Then the solution 
of Eqs. (HU) is 



(15) 



Q(t f )=Q Kb (t f ,t). 
Finally, using the function 

*/(*/.*) = H 2 (tf,t)^(t f ) = 2^ d TKb - 2 (r,t), (16) 

we can write the distribution in the form 

p = ^ 1 (tf,t)G(Q-Qf^ 1 (tf,ty i a f (tf,t)). (17) 



The function Q f k^ 1 (t / , t) has a simple meaning. Con- 
sider the noise-free trajectory that passes through the 
point Qf at time tf. This trajectory should have passed 
through the point QfK b ~ 1 (tf,t) at time t. As expected, 
the transition probability 1)17(1 is maximal for Q coincid- 
ing with this point. 

The function a 2 , Eq. I(l(j|) . is simply related to the 
function cr 2 (t), Eq. 1(10(1 . introduced earlier, 

*£(*) - i% 2 (tf,t)°b(tf) = °f(tf,t). (18) 
We will use this relation in what follows. 



B. General expression for the PPD near qb(t) 

From Eq. the PPD is determined by the product 
of the transition probability l(17|) and the ratio of the 
quasiperiodic distributions p(q,t)/ p(qf,tf). The distri- 
bution p(q,t) close to q b {t) was found earlier pun. it 

has the form 



Q = Q- Qb(t), 
Q F = u F /p b . (19) 



p(q,t) = / dpp'(p,Q,t), 
Jo 

£ f 1 
p'{p,Q,t) = -^=cxpj- — 

pKb(t, t') 



4>{p, t) = Qp\n 



[J-bh 



Here £ and t' are constants, s((f>) = s(<p + 2n) is a zero- 
mean 27r-periodic function, and Qf is the dimensionlcss 
modulation frequency. 

The function s((f>) in Eq. 1(19(1 plays the role of an in- 
stantaneous modulation-induced change of the activation 
energy. In the regime of strong synchronization of es- 
cape the minimal value s m of s(4>) satisfies the condition 
\s m \ D. The minima of s((f>) lie on the optimal es- 
cape paths [n| H3, p = pipt(i) = Po P t(t + tlt f ). Here 
n = 0, ±1, . . . enumerates periodically repeated MPEP's, 

see Fig. []]»; we set (gop t (*) = Popt(*)) = (<fo P t(»,.Popt (*))• 
Near the basin boundary q b {t) the optimal paths satisfy 
the linearized equations ((11|) and evolve in time as 

Popt(i) = «6" 1 (*,*')Popt(* / )j Qopt(i) = -Ob(t)Papt(t), 
Qopt(t) = goptCO " «.(*), IQoptl < Ag. (20) 

Eq. 120(1 describes how a given optimal path approaches 
(lb{t),p = 0) for t — ► oo. The parameter t' is determined 
by matching to one of the periodically repeated trajecto- 
ries 1(11(1 that start from (q a (t),p = 0) for t — > — oo. 

Expanding the function s(<f>) around its minima at 
the periodically repeated MPEP's, to second order in 
4>(p,t) — 0(popt(i),i) we obtain the probability distribu- 
tion 1(19(1 as a sum of contributions from the MPEP's, 

£e—~/ D ^ 



p'(p,Q,t) = 



D 



exp 



-r^(p,Q,t)/D 



,(«) 



pQ 



p 2 o 2 (t) , n; 



In 2 



(21) 
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Here s'^ = (d 2 s/d<f> 2 ) m is the curvature of the function 
s at the minimum. For strong synchronization, where 
\s m \ 3> D, we have sj[, > D as well, which is a conse- 
quence of s(</>) being a zero-mean periodic function. The 
quantity s m can be found [19l l2(| , along with the con- 
stant £ exp(— s m /D), by matching the periodic distribu- 
tion (|19fl to the distribution well inside the attraction 
basin. 

It follows from Eq. i|20|) that near the basin boundary 

P& k \t) =P$(t + kr F ) =M b - k p^(t). (22) 

With account taken of the relation Mb — exp(27r/f2F), 
this leads to the expression 



2irkn F s" In 



P 



27r 2 /c"s m - (23) 



Combining Eqs. JT7J and (|2~Tll and using Eq. (|T5j) we 

obtain for the PPD near the basin boundary, 

Ph(q,t\q f ,t f ) 

= ^ 1 (tf i t)G(Q-Qf^ 1 (tf,ty, < Tf{tf,t)) 
J2 n f QC dpexp[-r^(p,Q 1 t)/D] 



EnJ^dpexp [-rW(p,Q f ,t f )/D] 



(24) 



This expression gives the general form of the PPD near 
the basin boundary. We note that, even though the equa- 
tions of motion near qb(t) can be linearized, the PPD is 
generally non-Gaussian. This distortion is an important 
feature of escape dynamics. 



IV. ADIABATIC REGIME NEAR THE BASIN 
BOUNDARY 

We start the analysis with the case of slow modula- 
tion, Hp <C 1, where the motion can be described in the 
adiabatic approximation. As we show below, in the adi- 
abatic regime and for strong synchronization, 3> D, 
only one term contributes to each of the sums in Eq. (|24H . 
The shape of the PPD in this case is determined by the 
parameter 



n 2 F s'jD. 



(25) 



We call 9 the distortion parameter. This is because for 
9 <C 1 the escape current has a form of Gaussian peaks, 
whereas for la rger 9 the peaks of the current become non- 
Gaussian [T^ . Eol ]. Formally, 9 is equal to the ratio of 
two small parameters, the squared reduced modulation 
frequency ui f/ fib = &f and the noise intensity scaled by 
the effective modulation strength D/s'^. 

The physical meaning of 9 can be understood in the fol- 
lowing way. In the adiabatic picture one usually thinks of 
escape as occurring in the instantaneous potential U(q, t). 
Most likely it happens once per period at the time t m 
when the barrier height AU(t) = U(q b (t), t) - U(q a (t), t) 



is minimal, cf. Fig. ^ As explained in the Introduc- 
tion, the typical width of the time window for escape At 
is determined by the condition [d 2 AU / dt 2 ] m (At) 2 = D, 
where the subscript m indicates that the derivative is 
evaluated for t = t m . The parameter s m is (see Refl20() 

C = [d 2 AU/dt 2 ] m /Q 2 Ff i b {t m ) ~ [d 2 AU/dt 2 ] m /u; 2 F . 

Therefore the parameter 9 ~ [(At) 2 p, 2 ,] -1 is the squared 
ratio of the relaxation time of the system t T = p,^ 1 to At. 
It shows whether the system moves fast enough to escape 
while the barrier remains at its minimum or the barrier 
noticeably changes during escape leading to a delay of 
escape with respect to t m and a distortion of the tube of 
escape trajectories. 

Because the PPD evolves over time ~ t r , of interest 
is the time range \t — tf\ < t x <C Tp. Moreover, both 
t and tf should be close to t m . Then the instantaneous 
relaxation rate /Lt(,(t) can be approximated by its value 
fibm = Mb(^m)- The functions n b and a 2 become 

Kbm (*/>*) = cxp \pbm(tf - t)] , 

°%n (*/,*) = ff L {1 ~ exp [-2 M 6m(*/ - t)]} , (26) 



where 



a 2 b (t) « a 2 b (t m ) 



'bm 



i/m, 



6m- 



A. Weak distortion, 9 <C 1 

We consider first the limit of weak distortion, 9 <C 1 
(more precisely, the weak distortion condition has the 
form 01n 2 (s m /-D) <C 1, see below). In this limit one can 
think of motion in the fully adiabatic way, assuming that 
it occurs in a quasistatic potential U (q, t m ). The periodic 
distribution behind the attraction basin p(qf,t), which 
is proportional to the instantaneous escape rate, has a 
form of periodic in time Gaussian pulses, with width ~ 

(£>/O v v M- 

For small 9 expression (|24|l for the PPD can be simpli- 
fied. We will start with the analysis of the denominator 
in this expression. If there were no term ex Vips'^/ D = 9 
in r( n '(p,Qf,tf)/D, the typical values of p contributing 
to the integral over p would be ~ D/Qf. For such p there 
may be only one n for which the term 9\vi 2 [p/p^ t (tf)\ in 
j-W/D is small, whereas for all other n it is oc s'^/D [cf. 
Eq. 123|) ]. making the integrand exponentially small. A 
similar argument applies to the numerator in Eq. 12411 . 
except that, depending on Q, the typical values of p are 
of order of (D/^ bm ) 1 / 2 , D/Q, or -fi b mQ- 

Keeping only the leading term in the sums in the nu- 
merator and denominator of Eq. (|24|l and disregarding 
corrections oc 9, we obtain 

p h (q, t\q f ,t f ) =G(Q- QfK^(tf,t); cr/ m (*/,*)) 

Q 2 



erfc 



.(27) 
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In deriving this expression we also took into account that 
the final point qj is sufficiently far from the diffusion- 
dominated region behind the basin boundary, Qf ^> Id, 
and disregarded corrections ~ l D /Q"j. We note that the 
terms oc 8 in A 11 " 1 in the numerator and denominator in 
Eq. I|24[l have a logarithmic factor. This factor may be- 
come large in the weak-noise limit. One can show that 
for the corrections oc 8 to be small it is necessary that 
8 In 2 (sj^/ D) <C I. This condition is equivalent to the 
inequality At 3> ts, where At is the characteristic time 
window within which the barrier height is practically con- 
stant and ts is the Suzuki time, see Eq. I|28[l below. 

Equation (|27[1 for the PPD is further simplified for such 
times that the point QfK^ 1 (tf,t) is far behind the basin 
boundary compared to Ijj. This is the point that the 
noise- free trajectory arriving at Qf at time tf passes at 
time t. For such t, the PPD as a function of Q has 
the form of a Gaussian peak with variance DaJ(tf,t). 
This means that the trajectories arriving to the point 
Qf are close to the noise-free trajectory. The tube of 
these trajectories is diffusion-broadened, with width oc 
[D(t f - t)} 1 ' 2 for small t f ~ t. 

Because for 8 <C 1 escape as a whole occurs in the 
quasistatic potential U(q,t m ), the full calculation of the 
PPD described in Section VI leads to the same result as 
in the case of escape in a stationary potential studied 
earlier 15J. With increasing tf — t the peak of the PPD 
(|27[) crosses the basin boundary and enters the attraction 
basin. This is the slowest part of the PPD evolution. Its 
duration is determined by the Suzuki time 

t s = t I ln(C/D). (28) 

The peak is sharply broadened in this region. Deep on 
the intrawell side, — Q ^S> Id, the width of the peak be- 
comes independent of D, i.e. parametrically larger than 
the diffusion- limited width ~ Id- As tf — t increases fur- 
ther, the PPD peak approaches the attractor and narrows 
down, with the width becoming again diffusion-limited. 
This part of the evolution takes ~ t T . 



B. Strong distortion, 6 3> 1 

The shape of the PPD near the maximum changes dra- 
matically in the range where the modulation is still dy- 
namically slow, Clp <C I , but the distortion parameter 8 
becomes large, 8^1. Here, the shape of the potential 
barrier for escape changes as the particle crosses the dif- 
fusion region around the basin boundary. This leads to a 
strong change of the PPD compared to the picture based 
on the quasistatic barrier that was discussed before. 

For 8 ^ I in an important range oiQ,t and Qf,tf the 
integrands in the numerator and denominator in Eq. 1241) 

have sharp extrema as functions of p for p = p^ t (a more 
precise condition is specified below). Integration over p 
can be done by the steepest descent method. It gives the 
PPD in the form of a Gaussian distribution over Q with 



time-dependent center Q no and variance DaJ, 

p fl {q,t\q f ,t f ) = G{Q-Q m ;a f (t f ,t)), (29) 

where 

Qn = Q n (t\qf,tf) 

= Q$(*) + (Qf - K b\tf, t) (30) 

and no = no(g/,£/). Here we have taken into ac- 
count that the major contribution to each of the sums 
in Eq. (|24|1 comes from one term, with n — uq. 

It is clear from the analysis of the denominator in 
Eq. (|24|) that the value of no is determined by the condi- 
tion that p^l\tf)Qf/D ~ I for the noth optimal path. 
The very existence of such no follows from the fact that 
we consider tf for which the probability density p(qf,tf) 
is close to its maximum over tf, this maximum is reached 

once per period for p [ ^ t {tf)Qf/D = I 0[lj. The terms 
with |n — no | > I in the denominator in Eq. Q24f) are 

oc Ml l °-' n exp (-M™ "") < I 

where we have used Mf, = exp(27r/51i?) 3> I (this esti- 
mate is written for (t f)Q f / D = I). 

A simple qualitative argument shows that the same no 
gives a major contribution to the sum in the numerator 
of Eq. 124|) . Indeed, the terms with different n corre- 
spond to t changing by an integer number of the mod- 
ulation periods Tp, whereas in its central part (section 
C in Fig. ^) an optimal escape trajectory lasts for the 
time small compared to Tp. Therefore \t — tf \ <C Tp. 
The formal condition for this approximation to be true 
is that the Suzuki time ts <|28[) is small compared to the 
modulation period. 

The maximum of the Gaussian PPD peak H29|) lies at 
Q no (t\qf,tf). This function is a sum of the optimal path 
Qopti which is located inside the attraction basin, and 

the decaying in time term oc Qf — Q^pt(tf), which is de- 
termined by the final point (qf, tf) and is located outside 
the attraction basin. Its time dependence is particularly 
simple in the case where t f corresponds to the maximum 
of the distribution p(qf,tf), i.e., p^itf) — D/Qf, 

Q no (t\Q f ,t f ) « Qje-^t-*) ^- e M*»(t/-*). (31) 

The first term in Q na , Eq. I|3I|I . decreases with increas- 
ing tf — t, which describes approaching the basin bound- 
ary backward in time. The second term, on the other 
hand, increases with tf — t; this term describes the mo- 
tion, backward in time, from the boundary to the interior 
of the attraction basin. The motion at the boundary is 
initiated by noise, and therefore this term is oc D. For 
sufficiently long tf — t the distribution 

proaches the optimal path (t). The overall behavior 
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FIG. 3: (Color online) The reduced position of the maximum 
(left) and the reduced width (right) of the PPD in the 
adiabatic regime as a function of the reduced time fibm(tf —t). 



of Q no is shown in Fig. [21 It is seen that the motion is 
slowed down near the basin boundary, and the slowing 
down strongly depends on the noise intensity, in agree- 
ment with Eq. i|5T)l. 

The width of the Gaussian PPD peak l|29|) is diffusion- 
limited; its time dependence is given by Eq. (126[l and is 
also simple. It is shown in the right panel of Fig. [3J 

The condition that the integrands over p in Eq. 124|) 

are maximal for Pop" imposes a limitation on Q where 
the explicit expression for the PPD (|29[l applies. For Q 
close to the maximum of the distribution, \Q— Q no \ < Id, 

this condition takes the form 8 ^> IQopt (*)Ad|- There- 
fore expression l|29ll describes the distribution not only in 
the whole range between (Qf,tf) and the basin bound- 
ary, but also throughout the diffusion region around the 
basin boundary and a region that goes deeper into the 
attraction basin. 



V. NONADIABATIC REGIME NEAR THE 
BASIN BOUNDARY 

The PPD can display a qualitatively different behavior 
for nonadiabatic modulation, Qp ~ 1. With increasing 
tf — t the PPD inside the attraction basin can split into 
several peaks. Such splitting occurs already close to the 
basin boundary and is described by the general expres- 
sion J23J. Formally, the PPD has well-resolved multiple 
peaks when several terms in the numerator of Eq. i|24|) 
are of the same order of magnitude. 

In the nonadiabatic regime necessarily 8 1. There- 
fore integration over p in Eq. I|24[l can be done by the 
steepest descent method. The integrands are maximal 
for p — Po™].. Using Eq. ifTHl) one can show that each 
term in the numerator gives a Gaussian distribution over 
Q centered at Q n = Q n {t\qf,tf), Eq. (|3U[) . with variance 
Da 2 f (tf,t), Eq. I|lt)|) . This is similar to the adiabatic case 
(|29|l . except that now the PPD is a sum of appropriately 
weighted Gaussian peaks, 

n(at\a t\- En B n G (Q - Q n ; a f (t f ,t)) 
p h (q,t\q fl t f ) - . (32) 



The weighting factors B n are 

B n = B n (qf,tf) = x„exp(-x„), 

x n = x n (q f ,t f ) ^p ( ^ t {t f )Q f /D. (33) 

If we keep only one term in the numerator and denom- 
inator in Eq. (|32|l . with the same n = n$ given by the 
condition that B n is maximal for n — no, Eq. I|32|) goes 
over into Eq. (|29|l . However, in the nonadiabatic case 
there are regimes where the PPD as a function of {Q,t) 
may display several peaks. Eq. I|32[l gives their shapes 
near the maxima. 

The number of peaks of ph depends on Q,t, and ftp, 
as well as the final point (g/,i/). As before, we choose 
the final point so that the probability density behind the 
basin boundary p(qf,tf) is close to its maximum over 
tf, which occurs for x no {qf ,tf) ~ 1. The amplitudes of 
different peaks are 

Bn +k = x no M^ k exp (-x no M b ~ k ) , (34) 

with M;, = exp(27r/f2^) being the Floquet multiplier, 
Eq. (JJJ. In practice, in the whole range of modulation 
parameters where escape is strongly synchronized the fac- 
tor Mb is comparatively large. Therefore for x no ~ 1 the 
coefficients B ng +k rapidly decay with increasing |fc|, and 
only few peaks of the PPD can be simultaneously ob- 
served, primarily with non-negative k = 0, 1, . . .. 

The other condition for observing several peaks is that 
the distance between them exceed their width, that is in 
particular \Q no - Q no+ i| > Id- For Q ps Q no > l D , 
i.e., when the PPD maximum is well behind the basin 
boundary with respect to the attraction basin, we have 

\Q no -Qn +i\ - \Q ( Zt(t)\ - l 2 D^b(tf,t)/Q f « l D , where 
we have used that Q / 3> Id and that Kb(tf,t) ~ 1 when 
the PPD maximum is behind the boundary. Therefore 
in this region the PPD has indeed only one peak. 

For longer tf — t when Kb(tf,£) 1 and the peak po- 
sitions Qn : Qn +i are wei l inside the basin of attraction, 
the distance between them exceeds Id- In this range the 
corresponding peaks of ph can be resolved. The ratio 
of their amplitudes is given by the factor Mb, as seen 
from Eq. I|34|l . This is in good agreement with the re- 
sults of numerical simulations for a specific model shown 
in Fig. |3 as will be discussed below. 

We note that for flp 3> I synchronization of escape by 
modulation becomes weak. Although the system is far 
from the adiabatic limit, dynamics of the period-average 
coordinate is similar to dynamics in the case of a sta- 
tionary system. Even though the factor Mb becomes of 
order 1, of interest is the PPD with respect to the period- 
averaged coordinate, and this PPD does not display mul- 
tiple peaks. 

VI. THE PPD INSIDE THE ATTRACTION 
BASIN 

We assume throughout this Section that escape is 
strongly synchronized and 9 > 1. To find the PPD 
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for (q, t) deep inside the attraction basin we will use the 
equation 

Ph(q,t\q f ,tf) = J dq'p h (q,t\q',t')p h (q',t'\q f ,tf). (35) 

It applies for t < t' < tf, as follows from the 
definition J3J), and can be obtained from the 

Chapman-Kolmogorov equation for the transition prob- 
ability p{q f ,t f \q,t). 

As we show, t' in Eq. (|35[1 can be chosen in such a way 
that only a narrow range of q' contributes to the inte- 
gral. The corresponding q' are close but not too close to 
qb(t'). In this range both factors in the integrand, and 
hence the integral as a whole, can be explicitly calcu- 
lated. In particular, the function Ph{q' , t'\qj, tf) is given 
by Eq. 

The function ph(q,t\q' ,t') is the PPD inside the at- 
traction basin. For a periodically modulated system it 
was found earlier [2(|. Because in escape the system is 
likely to move close to one of the periodically repeated 
most probable escape paths q^,\(t), we are interested in 
Ph(q, t\q' , t') for both (q,t) and (q',f) lying close to such 
a path. If both (q, t) and (q' , t') are close to an nth path 

«7^™t(i), the corresponding PPD p^\q,t\q' ,t') is Gaussian 



pW(q,t\q',t') = G(q-q n (t\q',t');a n (t>,t)) (36) 

Here q n (t\q',t') is the value of the coordinate at time t 
on the optimal path that leads to the point (q' , t'). This 
optimal path is described by Eq. ifTTjl. 

The coordinate q n (t\q' , t 1 ) as a function of time is close 
to q^ t (t). We can seek it in the form q n (t\q',t') — 
9optW + 5q n {t\q' ,t'). To the lowest order in the devia- 
tion of (q' , t') from q^ t (t'), the function 6q n can be found 
from linearized equations (|llfl . This gives 

Sq n = V n (t)Sq n , 



where 



(37) 



V n {t) 



d qt K + -d qq {K 2 ) 



The boundary conditions for Sq n are 5q n (t'\q', t') = q' — 

9opt(^') an d limt^^oo 5q n {t\q' , t') = 0; the latter condi- 
tion simply means that the optimal path approaches the 
attractor q a (t) as t — > — oo. 

It follows from Eq. I|37|) that function f3 n (t) = 5q n /Sq n 
satisfies a first-order (Riccati) equation, 



(38) 



For q',q^\(t') both close to the basin boundary qb(t'), 
with account taken of the boundary conditions for Sq n 
the solution of Eq. I|37|l can be written in terms of f3 n as 

Sq n (t\q',t') = (Q 1 - Qll\(t')) Kn (t,t>), (39) 



K n (t,t') = exp 



dr/3 n (r) 



(40) 



We now discuss the asymptotic behavior of (3 n (t) . For 
not too large t' — t the function V n (t) in Eqs. l|57|l. 

should be calculated for q^ t {t) ~ qb(t)- For such V n (t), 
the solution 5q n (t) of Eq. l(3T|) can either exponentially 
increase or decrease with increasing t' — t. Of interest to 
us is the decreasing solution, which will ultimately go to 
zero for t' — t — > oo. A simple calculation shows that for 
this solution 



Pn(t) » fl h (t) = [d q K] qb(t) , 



{n) (t) 



Qopt 



q b (t). 



For much larger t' — t the optimal path <?optW a P~ 
proaches the attractor q a (t) , and then the function V n (t) 

in Eqs. l(3T|). iJSSJl should be calculated for q$ t (t) w q a (t). 
One can show that the solution 5q n (t) 
corresponds to 



for t 



Qa{t), 



-oo 



where cr 2 (t) is the variance of the periodic distribution 
about the attractor given by Eq. (|1UI) . 

To the lowest order in q' — q^\(t'), the variance of the 
prehistory probability distribution l|36l) is also expressed 
in terms of the function {3 n [2r| . 



dTK n 2 (T,t). 



(41) 



Using Eqs. ijlDjl. and (glj one can show that 

cr^(i',i — * — oo) = cr^(t). For t — ->■ — oo the distribu- 
tion Pfr l \q, t\q', t') goes over into the periodic Gaussian 
distribution p a (q,t), Eq. (JSJ), centered at the periodic at- 
tractor. 

It follows from Eqs. (JSJ-EHJ that, for q close to g£$(t), 

the PPD p^(g,t|g , ,t / ) displays a diffusion-broadened 
Gaussian peak as a function of q' with maximum close to 
9opt(^')- The displacement of the maximum over q' from 

9opt(*0 * s ^ fe - 9opt(*)]- On the other hand, the function 
Ph(q',t'\qf,tf) as given by Eq. is a sum of Gaus- 
sian distributions centered at the optimal escape paths. 
Therefore integration over q' in Eq. i|35[) can be done by 
the steepest descent method. The extreme values of q' 
are also close to the optimal escape paths. 

Different peaks of the PPD p h (q' ,t'\q f ,t f ), Eq. 
correspond to the MPEP's shifted by an integer num- 
ber of modulation periods. For each of these peaks one 
should use in Eq. the PPD p h (q, t\q't') given by 
with the appropriate n. 

The result of integration over q' describes the peaks of 
the PPD ph(q,t\qf,tf) for q close to (fcptW w hh appro- 
priate n. These peaks are Gaussian near the maxima, 
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Ph(q, t\q f ,t f ) = ^—^ — B n (q f ,t f )G(q - gg(t) - (q s - Q<g(f/)) K^C*/,*); a n (t f ,t)). (42) 

I 



In obtaining this equation we used the relations 

Kb(tf,t')K n (t',t) =K n (tf,t), 

a 2 n (t', t) + a 2 f (t f ,t')K- 2 (t', t) = a 2 n (t f ,t), 

which in turn follow from the relation (3 n (t) ~ fJ-b(t) in 
the harmonic region near the basin boundary. As a con- 
sequence, the result of integration over q' in Eq. (|35|l is 
independent of the matching time t' . 

Equation i|42[) is the central result of the present pa- 
per. It gives, in the explicit form, the PPD of finding the 
system at a position q inside the attraction basin at a mo- 
ment t, provided the system has been observed at a point 
qj outside the attraction basin at a time tf > t. It shows 
that, in the regime of strong synchronization, the peaks 
are Gaussian. They are centered at the most probable 
escape paths and are diffusion broadened throughout the 
attraction basin. This is in dramatic contrast with the 
PPD in the absence of synchronization, where the PPD 
peak inside the attraction basin is strongly asymmetric 
and its width is independent of the noise intensity D |l5j | . 

It follows from Eq. (|42[) that the PPD may have mul- 
tiple peaks inside the attraction basin. They can be ob- 
served only for strongly nonadiabatic modulation, where 
the Suzuki time is tg > Tp . In this case the system stays 
in the diffusion layer around the basin boundary qb{t) 
long enough to accumulate influxes from several period- 
ically repeated most probable escape paths q^\(t). On 
the other hand, the modulation period Tp should not 
be too short, because the strong synchronization of es- 
cape would be lost. Since synchronization loss occurs 
for t T ^> Tp, the PPD displays well resolved multiple 
peaks only in a limited parameter range. The limitation 
is more restrictive to the considered case where the final 
point (qf, tf) is close to the maximum of the distribution 
behind the barrier. If this condition is not imposed, the 
peaks are well resolved in a broader range, but measuring 
the PPD becomes more complicated on the whole. An 
example is discussed in the next Section. 

An important feature of the distribution (|42|l is weak 
dependence of the shape of the Gaussian peaks inside 
the attraction basin on the final point (qf,tf), which is 
a consequence of the smallness of the factor K _1 (i/,i). 
This shows that the PPD reveals the actual structure 
of the tubes of the paths followed in escape inside the 
attraction basin. In contrast, the relative amplitudes of 
the PPD peaks B n are sensitive to the choice of the point 



VII. RESULTS FOR A MODEL SYSTEM 

In this section, we present the results of numerical sim- 
ulations of the PPD for a simple model system and com- 
pare them with the analytical predictions. We consider a 
Brownian particle in a sinusoidally modulated potential 
of the form of a cubic parabola. The Langevin equation 
of motion has the form of Eq. with 

K{q,t) = q 2 - | +Acos(ujpt). (43) 

The dynamics was simulated using the second-order 
integration scheme for stochastic differential equations 
|27| . The system was initially prepared in the vicinity 
of the metastable state q a (t)- The final point (qf,tf) 
was chosen behind the basin boundary <#,(£). The PPD 
was calculated as a normalized probability distribution 
of paths q(t) arriving at the point qf for a particular 
modulation phase <pf — uptf (mod 2ir). 

An example of the full PPD is shown in Fig. The 
point (qf,tf) is chosen so that in escape the system is 
likely to pass near it (tf is determined modr^), that is 
the quasistationary distribution p(qf,tf) is close to its 
maximum over tf for a given qf behind the basin bound- 
ary [the parameter p^l 1 (tf)[qf — qb(tf)]/D is equal to 1.2, 
whereas the maximum of p(qf , tf) is expected where this 
parameter is equal to 1, 20]]. For the chosen modulation 
parameters and noise intensity the calculated activation 
energy of escape is R ss 0.0910 and R/D ss 9.1; the ra- 
tio of the modulation frequency to the relaxation rate 
ftp pa 2.24. We accumulated ~ 10 5 escape trajectories 
that arrive into a small area centered at (qf,tf), with 
width 5q = 0.02, St = 0.02r F . 

It is seen from Fig. that, in the regime of strong 
synchronization and for strongly non-Gaussian pulses of 
escape current, the peaks of the PPD are narrow both be- 
hind the basin boundary and inside the attraction basin. 
Moreover, for the chosen parameter values two distinct 
peaks of the PPD are well resolved inside the attraction 
basin. They correspond to the two paths the system is 
most likely to follow in escape. 

The positions of the PPD peaks on (q, t) plane are 
shown in Fig. 01 with full squares where there is one peak, 
and with crosses where two peaks are well resolved. They 
are compared with the periodically repeated optimal es- 
cape paths calculated by solving numerically the varia- 
tional equations (|llf> for the model , l|43|) . Such paths 
start from the periodic attractor q a (t) for t — > — oo and 
approach the basin boundary qb(t) for t — > oo, which are 
also shown in the figure. It is seen that, as expected 
from our analysis, the PPD maxima lie nearly on top of 
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FIG. 4: (Color online) The positions of the maxima of the 
PPD ph(q,t\qf ,tf) in Fig. 01 which show the most probable 
paths followed by the system in escape. The data of sim- 
ulations are shown by full squares where the PPD has one 
peak and by crosses where two peaks are well resolved. Solid 
lines show periodically repeated most probable escape paths 
?opt(*) f° r the model (0, (1431) . Dashed lines show the basin 
boundary qt(t) and the attractor q a (t). 



the most probable escape paths. The small deviation is 
due to diffusion broadening and associated small asym- 
metry of the PPD peaks for the noise intensity used in 
the simulations. Fig. |I] demonstrates that studying the 
PPD provides a direct way of observing most probable 
escape paths. 

The observed shape of the PPD peaks is compared with 
the theory in Fig. El which shows the cross-sections of the 
PPD at several instants of time counted off from the final 
time tf. For small tf — t the system is behind the basin 
boundary and moves close to the noise-free trajectory 
leading to (g/,t/) . The top left panel refers to the case 
where the system is localized close to the boundary. Here 
and for smaller tf — t the PPD has a single sharp peak. 

For earlier time (larger tt — t) the system could ei- 
ther be moving towards the basin boundary along the 
most probable escape path or could have been fluctuat- 
ing about the basin boundary after it had arrived to its 
vicinity along the previous MPEP (shifted by rp). As 
we showed analytically, the probability of staying near 
the boundary is smaller, but the PPD may still display 
two peaks. This is seen in the right top panel. The main 
peak corresponds to motion along the MPEP, which is 
close to qb(t) for the chosen time. The higher-q shoulder 
corresponds to the poorly resolved (for the chosen time) 
peak for fluctuations about the basin boundary. 

For still larger but not too large tj — t the escaping sys- 
tem should have been moving towards the basin bound- 
ary. In the present case it most likely followed one of the 
two periodically repeated MPEPs, with different proba- 
bilities. Well-resolved peaks of the PPD in this range are 
seen in the left lower panel of Fig.JSJ For large tf — t com- 
pared to the relaxation time and rp, the system should 
have been fluctuating about the attractor q a (t). The 




FIG. 5: (Color online) Cross-sections of the PPD 
Ph(q,t\qf ,tf) shown in Fig. [5] as functions of the coordinate 
q for the time tf — t = 0.37>, 0.5l>, t>, 27>. The results re- 
fer to the model system (Jy, l|43|l with A = 0.7, ujf = 2, 
D = 0.01, qf = 0.8, tf = 0.5tf- The point (qf,tf) is close to 
the expected maximum of the distribution behind the basin 
boundary. Solid lines show the expression 1421 . Squares show 
the results of simulations. 



PPD in this case should have a single peak, which is 
seen in the right lower panel. 

It is seen from Fig. [5] that the results of simulations 
agree with the analytical results. Not surprisingly, the 
observed peaks are broader than the asymptotic theory 
predicts. This is a consequence of the relatively large 
noise intensity used in the simulations. 



VIII. CONCLUSIONS 

In this paper we have studied the prehistory prob- 
ability distribution for activated escape in periodically 
modulated systems. We have shown that the PPD as a 
function of coordinates can display one or several narrow 
peaks within the basin of attraction to the metastable 
state. They correspond to narrow ridges of the PPD 
in (q,t) space. The ridges are centered at the periodi- 
cally repeated most probable escape paths. Their cross- 
sections are Gaussian near the maxima and are diffusion 
broadened, see Eq. (|42p. 

The number of the PPD peaks that can be observed de- 
pends on the parameters of the system. The multi-peak 
structure is best resolved in a limited parameter range. 
On the one hand, the modulation should not be too slow, 
so that the system does not follow it adiabatically. On 
the other hand, it should not be too fast, so that escape 
events are strongly synchronized and the system dynam- 
ics is not described by period-averaged coordinates. 

Most of the results of the paper refer to the case where 
the final point {qf,tf) is chosen so that the escaped sys- 



12 



tern has a high probability density of passing through 
this point. Such choice simplifies the experimental ob- 
servation of the PPD. For the corresponding (q/,tf) the 
amplitudes of different peaks of the PPD [heights of the 
ridges in (q, t) space] differ from each other significantly, 
as seen from Eq. ()34|l and Figs. El El These amplitudes 
are sensitive to the choice of (qf, tf). In contrast, the po- 
sitions and shapes of the PPD peaks very weakly depend 
on (q/,tf). This shows that the PPD provides a means 
for studying the distribution of trajectories leading to 
escape. 

We have performed extensive numerical simulations of 
escape for a model Brownian particle in a modulated po- 
tential well. The simulations confirm the possibility to 
clearly observe most probable escape paths. They are in 
a good qualitative and quantitative agreement with the 
analytical theory. 

Observing most probable escape paths is not only in- 
teresting on its own, but also has broader implications. 
First, periodically modulated systems are an important 
class of systems far from thermal equilibrium. In contrast 
to the case of equilibrium systems, optimal fluctuational 
paths in nonequilibrium systems have no immediate re- 



lation to dynamical trajectories in the absence of noise; 
in particular, they may not be obtained by just reversing 
time. The MPEP's can display interesting and counterin- 
tuitive behavior 0, [2^ , and studying them provides an 
insight into general features of dynamics away from ther- 
mal equilibrium. Second, understanding the dynamics of 
a system in escape paves the way for efficient control of 
this dynamics and the escape probability itself; as will be 
discussed in a separate publication, such control can be 
accomplished by comparatively weak field pulses applied 
in the right place and at the right time. 

In conclusion, the results of this paper suggest a way 
of direct observation of most probable escape paths, in 
space and time. They also describe, qualitatively and 
quantitatively, the distribution of the trajectories fol- 
lowed in escape. 
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